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Abstract 

?*■ Accurate detection of mitosis plays a critical role in breast cancer histopathology. 

w Manual detection and counting of mitosis is tedious and subject to considerable 

^ inter- and intra-reader variations. Multispectral imaging is a recent medical imag- 

'— ' ing technology, proven successful in increasing the segmentation accuracy in other 

T_H fields. This study aims at improving the accuracy of mitosis detection by develop- 

^ ing a specific solution using multispectral and multifocal imaging of breast cancer 

.^ histopathological data. We propose to enable clinical routine-compliant quality 

O of mitosis discrimination from other objects. The proposed framework includes 

comprehensive analysis of spectral bands and z-stack focus planes, detection of 
expected mitotic regions (candidates) in selected focus planes and spectral bands, 

cF) computation of multispectral spatial features for each candidate, selection of mul- 

tispectral spatial features and a study of different state-of-the-art classification 
methods for candidates classification as mitotic or non mitotic figures. This frame- 

^ work has been evaluated on MITOS multispectral medical dataset and achieved 

60% detection rate and 57% F-Measure. Our results indicate that multispectral 
spatial features have more information for mitosis classification in comparison 
with white spectral band features, being therefore a very promising exploration 
area to improve the quality of the diagnosis assistance in histopathology. 
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1. Introduction 

Breast Cancer (BC) is the most commonly diagnosed cancer after the skin 
cancers, and the second leading cause of cancer death, following lung cancer 
among U.S. women [25]. In 2012, an estimated 226,870 new cases of invasive 
BC and 39,510 BC deaths in women are reported in U.S. In addition, BC inci- 
dence and mortality rates have been increasing rapidly in economically less de- 
veloped countries. According to World Health Organization, the reference process 
for breast cancer prognosis is histologic grading that combine tubule formation, 
nuclei atypia and mitotic counts [|2l|9|. This assessment of tissue sample is syn- 
thesized into a diagnosis that would help the chnician determine the best course of 
therapy. Several CAD solutions exist for the detection of tubule formation Il30ll29ll 
and nuclei atypia [|6l El 111 ISl but only few are dedicated to mitotic counts (MC) 

Ellll. 

In histopathology, H&E is a well-established staining technique that exploits 
intensity of stains in the tissue images to quantify the nuclei and other structures 
related to cancer developments [HI. In this context, image -processing techniques 
are devoted to accurate and objective quantification and localization of cancer evo- 
lution in specific regions of the tissue such as cytoplasm, membranes and nuclei 
[|28ll . From the chromatic viewpoint, nuclear regions are characterized by non- 
uniform stain intensity and color, thus preventing a trivial classification based on 
color separation. In addition, the superposition of tissue layers, as well as the 
diffusion of the dyes on the tissue surface, may bring the stains to contaminate 
the background or other cellular regions, which are different from their specific 
target. 

One of the most difficult fields in histopathological dataset analysis is spa- 
tial analysis, more specifically automated nuclei detection and classification [[TT|. 
The objective of nuclei classification is to assign different labels to different type 
of nuclei as normal, cancer, mitotic, apoptosis, lymphocytes etc. In addition, 
quantitative characterization is important not only for clinical applications (e.g., 
to reduce/eliminate inter- and intra-observer variation in diagnosis) but also for re- 
search applications (e.g., to understand the biological mechanisms of the disease 
process) iTlSl . 

Image analysis in cytology has been studied for years and numerous solutions 
[|36l [34l m im [T4|| have thus been proposed in the literature. The application 
of these solutions to histopathology is rather complicated due to the radical dif- 



ferences between the two imaging modalities and to the highly complex image 
characteristics. Indeed, in the case of histology images, cellular structures and 
functions are studied embedded in the whole tissue structure, presenting various 
cells architecture (gland formation, DCIS) very difficult to handle with usual pat- 
tern recognition techniques. Nevertheless, recent work [|26l[T7l[T8l show great po- 
tential for Computer Assisted Diagnostic of histopathological datasets for breast 
cancer diagnosis. 

Multi spectral imaging (MSI) has the advantage to retrieve spectrally resolved 
information of a tissue image scene at specific frequencies across the electro- 
magnetic spectrum. MSI captures images with accurate spectral content corre- 
lated with spatial information and reveals the chemical and anatomic features of 
histopathology [|20li22| . This modality provides option to biologists and patholo- 
gists to see beyond the RGB image planes that they are accustomed to. Recent 
publications ifTOl l23l l37l \T9 ] have begun to explore the use of extra informa- 
tion contained in such spectral data. Specifically, there have been comparisons 
of spectral methodologies which demonstrate the advantage of multispectral data 
[|2Tl[T3l . The added benefit of MSI for analysis in routine H&E histopathology, 
however, is still largely unknown, although some promising results are presented 
in [l32l[T0l[T9ll37ll . As far as we know, there is no existing study of the advantage, 
or lack of thereof, of MSI for automation of MCs in breast cancer histopathology. 
We propose here to extend the already successful work of [[IB] to support MSI and 
illustrate the advantages, paving the way for a better, automated detection of MCs 
in BC histopathology. 

The reminder of the paper is organized as follows. Section |2] reviews the 
state-of-the-art multispectral methods, particularly in object or region detection 
in histopathology, related to this research work. Section [3] describes the proposed 
framework for mitotic figure detection. Experiment and results are presented in 
section]?} Section [5] contain the discussion part. Finally, the concluding remarks 
with future work are given in section |6j 

2. Literature Review 

The main idea for extracting texture features from MSI is the use of combined 
spectral and spatial information for discrimination of regions or objects. We found 
few methods in the MSI literature for texture characterization of histopathological 
images. Some of them employed single band of MSI ll27l [37l l26l and other used 
multiple bands of MSI [|19[ [3||. Fernandez et al. [jlOl coupled high-throughput 
Fourier transform infra-red (FTIR) spectroscopic imaging of tissue microarrays 



with statistical pattern recognition of spectra indicative of endogenous molecular 
composition and demonstrate histopathological characterization of prostatic tis- 
sue. They explicitly defined metrics consisting of spectral features that have a 
physical significance related to tissue biochemistry and facilitating the measure- 
ment of cell types. 

Recently, Masood et al. [27] proposed a colon biopsy classification method 
based on spatial analysis of hyper- spectral image data from colon biopsy samples. 
Initially, using circular local binary pattern algorithm, spatial analysis of patterns 
is represented by a feature vector in selected spectral band. Later, classification 
is achieved using subspace projection methods like principal component analy- 
sis, linear component analysis and support vector machine. Recently, Wu et al. 
[|37ll proposed a multilayer conditional random field model using a combination 
of low-level cues and high-level contextual information for nuclei separation in 
high dimensional data set obtained through spectral microscopy. In this approach, 
the multilayer contextual information was extracted by an unsupervised topic dis- 
covery process from spectral images of microscopic specimen, which efficiently 
helps to suppress segmentation errors caused by intensity inhomogeneity and vari- 
able chromatin texture. Malon et al. [l26l demonstrated a segmentation based 
features with convolutional neural networks using selected focus plane and spec- 
tral band for identification of mitotic figures in breast cancer histopathology and 
achieved best classification accuracy (F-measure = 59%) on multispectral dataset 
during ICPR contest 2012 [|33l. These approaches [|27l|371|26| are limited to sin- 
gle spectral band by ignoring additional relevant information from other spectral 
bands. 

Boucheron et al. [[3| presented an analysis that the additional spectral bands 
contain further information useful for nuclear classification in histopathology as 
compared to the three standard bands of RGB imagery. Using all image bands, 
they reported 0.79% improved in performance as compared to next best perform- 
ing image type. This method is limited to pixel level features. One possibility 
of additional improvement in object classification is multispectral spatial analysis 
using texture features. 

Khehfi et al. [|T9l proposed a spatial and spectral gray level dependence 
method in order to extend the concept of gray level co-occurrence matrix by as- 
suming the presence of texture joint information between spectral bands. Having 
the fact that some spectral bands have more relevant information for specific ob- 
ject or region classification as compared to other spectral bands. This approach 
is limited to single spatial feature as computed from all spectral bands rather than 
one different spatial feature for each spectral band. 



Each object or region has different level of relevant information. In the pro- 
posed methodology, we address the shortcomings of previous works, including 
(1) comprehensive analysis of multispectral spatial features (MSSF) in all bands 
rather than single band [271 [571 |26l and (2) extracting MSSF in order to discrim- 
inate mitotic figures from other nuclei and microscopic objects. The main novel 
contributions of proposed work are (1) selection of z-stack focus planes and mul- 
tispectral bands for detection of candidates; (2) a multispectral spatial and mor- 
phology features computation which leverages discriminant information from a 
given nucleus across bands and (3) a selection of features from multispectral fea- 
tures vector for classification of mitotic figures in breast cancer histopathological 
images. 

3. Proposed Framework 

3.1. Dataset 

We evaluated the proposed framework on multispectral MITOS dataset OSlI 
a freely available medical mitosis dataset. The data set is made up of 50 high 
power fields (HPF) coming from five different slides scanned at 40X magnification 
using a 10 bands multispectral microscope. There are 10 HPFs per slide and each 
HPF has a size of 512 x 512/im^ (that is an area of 0.262 mm^). The spectral 
bands are all in the visible spectrum. In addition, for each spectral band, the 
digitization has been performed at 17 different focus planes (17 layers Z-stack), 
each plane being separated from the other by 500 nm. For one HPF, there are 
170 gray scale images (10 spectral bands and 17 layers Z-stack for each spectral 
band). These 50 HPFs contain a total 322 mitotic figures. The training data set 
consists of 35 HPFs containing 224 mitotic figures and evaluation data set consists 
of 15 HPFs containing 98 mitotic figures [|33l . Fig. [T] shows the spectral coverage 
of each of the 10 spectral bands of the multispectral microscope. 

3.2. Proposed Framework 

In this paper, we propose a framework for MC in Multispectral BC histopathol- 
ogy as shown in Fig. [2j A comprehensive analysis of all multispectral bands 
(10 bands) and z-stack focus planes is performed for detection of mitosis candi- 
dates. Then candidates are detected using thresholding and morphological oper- 
ations on selected band and focus plane. A multispectral features vector is com- 
puted for detected candidates having intensity and texture features across all bands 
of multispectral images. In addition, using segmented regions of detected candi- 
dates, morphological features are also computed. A feature selection algorithm 
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Figure 1 : Spectral bands of the multispectral microscope and examples for each 
band. 



is employed on this multispectral features vector in order to save the computation 
cost, to discard any redundancy in the data, and to improve classification accuracy. 
Classification is achieved using Bayesian, decision tree (DT), neural network as 
well as linear and non-hnear support vector machine (SVM) classifiers. A side 
advantage of performing the spatial analysis on a multiple band simultaneously 
is to investigate whether improvement in accuracy can be achieved with carefully 
selected multispectral features over those methods [.271 l37l |26J which use single 
band data. 

3.2.1. Candidate Detection 

Initially, we perform gradient analysis on all the focus planes. Planes with 
lower average gradient (i.e. out of focus) are discarded. We select focus planes 5, 
6, 7, 8, 9 and 10 (total 6 focus planes). On selected focus planes, we compute the 
spectral responses of mitotic nuclei and background regions for all the available 10 
spectral bands as represented in Fig. [3] while responses of mitotic and non mitotic 
nuclei are presented in Fig.|4j Note that spectral band 1 (white band), in nature, is 







j<- 




















Nw 




Candidate Detection 








,. J, 




f \ 
























1 


-*■ 


Focus Plane 
Selection 




Spectral Band 
Selection 




Thresholding 
& Morphology 




Candidate 
Selection 


1 














r 


















Marker 
Image ^ 


— Mask 


, 












"~\ 






> 












Spatial Features 




Morphological 
Features 






Feature 
Selection 




Classification 












10 bands, 17 layers Z-stack 
Multlspectral Data 


Multispeetral Features Computation 


Features Selection and Classification J 



Figure 2: Proposed Framework. 

different from others and might serve as reference as it contains all the information 
that other bands are containing, even if at a lower resolution. The peaks of the 
mitotic and non-mitotic nuclei are almost similar. Mitotic nuclei have different 
peaks than background regions. Multispeetral data is thus able to differentiate 
between different tissue parts but spectral response of mitotic and non-mitotic 
nuclei is not distinguishable. The process of cell division has four different stages 
and each has different shape, size and textures. This motivates further spatial and 
morphological analysis on multispeetral data to achieve reasonable classification 
of regions into mitotic and non-mitotic types. 

Using selected focus plane and spectral band, we eventually perform thresh- 
olding followed by morphological processing to eliminate small regions and fill 
holes and later select the candidates by filtering based on minimum size (37/um^ 
i.e., 200 pixels) and maximum size (lOOOyum^ i.e., 5405 pixels) of mitotic nuclei 
while keeping the number of candidates to classify lower. 

S.2.2. Multispeetral Features Computation 

Instead of single band as in [[3] |27l |37l |26'| we compute MSSF vector having 
intensity and textural features in all bands (10 bands). In addition, we also include 
morphological features (such as area, roundness, elongation, perimeter and equiv- 
alent spherical perimeter) that are computed from regions segmented during can- 
didate detection. The morphological features reflect the phenotype information of 
mitotic nuclei. Using spatial information in multispeetral bands, we compute five 
intensity-based features including mean, median, variance, kurtosis and skewness 
for candidate. This results in 50 multispeetral intensity-based features. Haralick 
co-occurrence (HC) lfT6l and run-length (RL) [12J features are computed with one 
displacement vector in four direction (0", 45°, 90°, 135°) for all the spectral bands 
as in [il8il . These multispeetral textural features are rotationally invariant. So by 
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Figure 3: Histogram analysis of mitotic and background regions in 10 spectral 
bands. 
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Figure 4: Histogram analysis of mitosis and non- mitosis regions in 10 spectral 
bands. 



making average in all four directions, eight HC and ten RL features are computed 
for each candidate in each single spectral band. The resulted multispectral textural 
features vector consists of 80 HC features and 100 RL features for each candidate. 
The final MSSF vector contains 235 features for each candidate. 

3.2.3. Feature Selection and Classification 

Conceptually, a large number of descriptive features are highly desirable for 
classification of nuclei as mitosis or non-mitosis types. However, when using all 
computed features (i.e. 235 features) the classification performance is poor. Some 
features are irrelevant for classification and some features are redundant degrading 
the classification performance. The consistency subset evaluation method [24J is 
employed to select a subset of features that maximize the consistency in the class 
values. A projection of subset of features is performed on training dataset in order 
to evaluate the worth of subsets of features by the level of consistency in the class 
value. The consistencies of these subsets are not less than that of the full set of 
features. Later, we used these subsets in conjunction with a hill climbing search 
method, augmented with backtracking value 5, which looks for the smallest subset 
with consistency equal to that of the full set of features. This procedure achieved 
88% reduction in the dimensionality of features set by selecting 28 features. 

The selected features contained one morphological, 7 intensity and 20 textural 
features in different spectral bands except spectral band 5. The selected features 
set is used to train different classifiers like decision tree (Funcational Tre-FT), 
neural network (Multilayer Perceptron-MP), Bayesian, linear SVM (L-SVM) and 
non-linear SVM (NL-SVM) classifiers [35]. Throughout the experiments, the 
parameters used in FT classifier are numBoostinglterations =15 and minNumln- 
stance = 15, in MP are learning rate = 0.3 and momentum = 0.2, in L-SVM classi- 



fier are bias = 1, cost = 1, eps = 0.01 and kernel = L2-loss SVM and in NL-SVM 
classifier are kernel = rbf, degree of kernel = 3, eps = 0.001 and loss = 0.1. 

4. Experiment Results 

The proposed framework is evaluated on MITOS multispectral dataset [|38l . 
The results of candidate detection and classification methods are compared with 
ground-truth information provided along with the dataset. The metrics used to 
evaluate the mitosis detection of each method include: number of true positive 
(TP), number of false positives (FP), number of false negative (FN), sensitivity 
or true positive rate (TPR), precision or positive predictive value (PPV) and F- 
Measure. In addition to MITOS contest metrics, the proposed framework is also 
evaluated with other state of the art method as 5-fold cross validation. 

4.1. Candidate Detection in Different Spectral Bands 

To gain a better understanding of the relative contributions of specific spec- 
tral bands, we perform candidate detection in all available spectral bands and six 
selected z-stack focus planes. The results of candidate detection are ranked ac- 
cording to F-Measure and reported top three rank results in Fig.[5](a). The focus 
planes six and seven have more information for candidate detection as compared 
to other focus planes and a bar graph representing candidate detection in focus 
planes six and seven with respect to F-Measure is shown in Fig.|5](b). 

The important parameter that affects the classification is the unbalance training 
set having a large number of non-mitosis as compared to few mitosis. According 
to highest F-Measure in all spectral bands and focus planes, the number of de- 
tected candidates in training and evaluation sets is shown in Fig. [5](c). Spectral 
band eight in focus plane six (590-640 nm) is selected for final candidate detection 
as it detect less FP as compared with other spectral bands and more TP than other 
spectral bands as well. On training and evaluation sets, the candidate detection 
using spectral band 8 in focus plane 6 detects 3583 and 1655 candidates, contain- 
ing 203 and 91 ground-truth (GT) mitosis from a total 224 and 98 GT mitosis, 
respectively. Therefore, among the entire detected candidates, there are 3380 and 
1564 non-mitosis in the training and evaluation sets, respectively. The candidate 
detection stage generates a large number of non-mitosis and misses 21 and 7 GT 
mitosis, from training and evaluation sets respectively. 
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(a) Top three rank focal planes. 
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Table 1: Classification result (GT = 98) on evaluation set. 



Experiments 


TP 


FP 


TPR 


PPV 


F-Measure 




Using 
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L-SVM Classifier 


33 
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NL-SVM Classifier 


40 


28 
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48% 


Bayesian Classifier 
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4.2. Candidate Classification 

4.2.1. Experiment 1: Classification using separate training and evaluation sets 
In experiment 1, training set is used to train the five selected classifiers and 

evaluation set is used to test the classification accuracy of proposed framework 
as shown in Table [TJ When all MSSF are used for training and evaluation, we 
get good TPR and F-Measure with Bayesian and FT classifiers. The MP, L- 
SVM and NL-SVM classifiers have few FP but also few TP and resulted low 
F-Measure. In case of selected MSSF, F-Measure is increase more in Bayesian 
classifiers with highest TPR than other classifiers. Overall, Bayesian classifier 
reports highest TPR (60%) and F-Measure (57%) while FT classifier has second 
highest F-Measure with few FP and high PPV. 

4.2.2. Experiment 2: Classification using 5 -fold cross validation 

In experiment 2, the assessment of classification performance is made using 
5-fold cross validation by combining both training and evaluation set as shown 
in Table [2} In case of full set of MSSF during training and evaluation, FT clas- 
sifier outperforms with highest TPR and F-Measure. Overall, Bayesian classifier 
reported highest TPR (62%) but more FP as well. In case of selected MSSF, as 
compared with L-SVM, NL-SVM and MP classifiers, FT and Bayesian classifiers 
report better F-Measure 54% and 52%, respectively. 
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Table 2: Classification result (GT = 322) using 5-Fold Cross Validation. 



Experiments 


TP 


FP 


TPR 


PPV 


F-Measure 




Using 


aU MSSF 






L-SVM Classifier 


95 


30 


30% 


76% 


43% 


NL-SVM Classifier 


95 


37 


30% 


72% 


42% 


Bayesian Classifier 


201 


384 


62% 


34% 


44% 


MP Classifier 


132 


90 


41% 


59% 


49% 


FT Classifier 


168 


140 


52% 


55% 


53% 




Using selected MSSF 






L-SVM Classifier 


96 


24 


30% 


80% 


43% 


NL-SVM Classifier 


99 


35 


31% 


74% 


43% 


Bayesian Classifier 


169 


165 


52% 


51% 


52% 


MP Classifier 


136 


95 


42% 


59% 


49% 


FT Classifier 


141 


62 


44% 


69% 


54% 



4.2.3. Experiment 3: Classification using white spectral bands vs other multi- 
spectral bands 

To investigate the relative contribution of multispectral bands, we also perform 
a comparative study of mitosis classification using spatial features from white 
spectral band (spectral band 1) with spatial features from multispectral bands and 
the achieved results are shown in Table |3] It is important to note that the multi- 
spectral bands features (MSBF) excluding white spectral band outperforms with 
all classifiers and reports highest F-measure 55% with FT classifier. The classifi- 
cation results obtained with all five classifiers are worst using white spectral band 
features (WSBF). This experiment illustrates that multispectral band have much 
more information for mitosis classification than white bands. 

4.2.4. Experiment 4: Multispectral texture features vs multispectral intensity fea- 
tures 

In this experiment, we explore the impact of multispectral texture features 
(MSTF) and multispectral intensity features (MSIF) on the mitosis classification. 
The results are shown in Table |4j In case of MSIF, we get high F-Measure 55% 
using MP classifier. The highest F-Measure (56%) is achieved with MSTF with 
FT classifier. These results illustrate that MSTF have more information for mitosis 
classification as compared MSIT. 
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Table 3: Classification result (GT = 322) using 5-Fold Cross Validation. 



Experiments 


TP 


FP 


TPR 


PPV 


F-Measure 




Usin 


gWSBF 






L-SVM Classifier 


25 


13 


8% 


66% 


14% 


NL-SVM Classifier 


28 


27 


9% 


51% 


15% 


Bayesian Classifier 


84 


236 


26% 


26% 


26% 


MP Classifier 


35 


40 


11% 


47% 


18% 


FT Classifier 


35 


36 


11% 


49% 


18% 




Usin 


gMSBF 






L-SVM Classifier 


82 


57 


25% 


59% 


36% 


NL-SVM Classifier 


83 


78 


26% 


52% 


34% 


Bayesian Classifier 


197 


615 


61% 


24% 


35% 


MP Classifier 


109 


43 


34% 


72% 


46% 


FT Classifier 


140 


84 


43% 


63% 


51% 



Table 4: Classification result (GT = 322) using 5-Fold Cross Validation. 



Experiments 


TP 


FP 


TPR 


PPV 


F-Measure 




Usin 


ig MSIF 






L-SVM Classifier 


81 


38 


25% 


68% 


37% 


NL-SVM Classifier 


86 


63 


27% 


58% 


37% 


Bayesian Classifier 


154 


278 


48% 


36% 


41% 


MP Classifier 


156 


91 


48% 


63% 


55% 


FT Classifier 


140 


83 


43% 


63% 


51% 




Usin 


gMSTF 






L-SVM Classifier 


84 


41 


26% 


67% 


38% 


NL-SVM Classifier 


88 


69 


27% 


56% 


37% 


Bayesian Classifier 


200 


553 


62% 


27% 


37% 


MP Classifier 


137 


65 


43% 


68% 


52% 


FT Classifier 


167 


105 


52% 


61% 


56 %c 
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5. Discussion 

The results seem to indicate the best scores are achieved using the same focus 
plane, or couple of focus planes, across bands. Each of those planes is exhibiting 
maximum gradient intensity above the rest of the focus planes. It means that in 
further experiments we could locate the best focus planes first, and then apply our 
framework on that limited dataset. In other words, finding the best focus plane 
and finding the best bands are separable problems. 

The best F-Measure for candidate detection is achieved on spectral bands 8, 
9 and 1 respectively. The fact that candidate detection using bands 8 and 9, the 
proposed framework achieves better results than when using the full spectrum 
(band 1) supports the claim that MSI improves the accuracy of the framework. As 
bands 8 and 9 actually overlap in terms of spectrum, it would be interesting to try 
to apply spectral unmixing between bands 8 and 9 to see if it can further improve 
the results. The results illustrate clearly the improved accuracy resulting of the 
selection process. 

We use the analysis of different subset of multispectral features as complement 
to the analysis of performance on MSSF. Specifically, 5-fold cross validation is 
employed to perform comprehensive analysis of proposed framework. In order to 
study the MSSF for classification of mitosis figures, we perform features selec- 
tion in different spectral bands by studying which spectral bands have redundant 
and irrelevant information for mitosis classification. Spectral band 5 is irrelevant 
for mitosis discrimination. Most of selected MSSF belong to spectral band 7 and 
8. Most of the selected features are texture in different spectral bands, especially 
spectral bands 7 and 8, which strengthening the importance of multispectral spa- 
tial analysis. 

Specifically we have shown in experiment 3 that most of the information is 
present in multispectral band excluding white spectral band that is not suitable for 
mitosis classification. According to experiments 4, multispectral texture features 
are more helpful for mitosis discrimination as compared to multispectral intensity 
features. 

Separate training and evaluation sets as provided with MITOS dataset ll38ll 
have been used for training and evaluation of proposed framework. To put the 
quality of this result in perspective, in comparison with MITOS contest result [38| 
proposed method with selected MSSF and Bayesian classifier managed to achieve 
second highest TPR and F-Measure. The comparison of proposed framework 
results with MITOS contest results are shown in Fig. [6} Fig. |7] illustrates the ROC 
curve obtained with selected MSSF with FT classifier. This clearly demonstrates 
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Figure 6: Comparison of proposed framework results with MITOS contest result. 

that our new proposed framework results in an improved ability to distinguish 
mitosis from other objects. The proposed framework should be considered as the 
state of the art. 

6. Conclusion 

An automated mitosis detection framework for breast cancer MSI based on 
multispectral spatial features has been proposed. Initially, candidate detection 
is performed in selected spectral band and z-stack focus plane. Then, we com- 
pute multispectral spatial features for each candidate, a highly efficient model for 
capturing texture features for region (nuclei) discrimination. The proposed frame- 
work reaches the best levels of the MITOS contest results. Both features set (MSIF 
& MSTF) shown similar classification performance for mitosis. In future work, 
we plan to investigate unmixing of bands and/or better selection of (sets of) bands 
of interest. The pre-selection of the focus plane (or volumes) is also of great 
importance to reduce the complexity of the dataset and improve the actual per- 
formances to reach clinical operational acceptance expected by our professional 
consortia. 
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Figure 7: The ROC curve using selected MSSF with FT classifier. 
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